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, Abstract. In a lattice formulation of 5 f [/(2)-gmodynamics, the spontaneous 

O'N \ generation of chromomagnetic fields at high temperature is investigated. A procedure 

to determine this phenomenon is developed. By means of the % 2 -analysis of the data 
set accumulating 5 — 10 x 10 6 Monte Carlo configurations, the spontaneous creation of 
d ' the Abelian color magnetic field is indicated. The common generation of the magnetic 

field and Ao-condensate is also studied. It is discovered that the field configuration 
Oh, consisting of the magnetized vacuum and the Ao-condensate is stable. 
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1. Introduction 



t^. \ Nowadays it is generally accepted that the non-linearity of non-Abelian gauge fields 
could result in formation of field condensates. The first model of gluon condensate, so 
called " color ferromagnetic vacuum state", was proposed thirty years ago by Savvidy pQ. 
It describes the spontaneous generation of the uniform Abelian chromomagnetic field 
H = const due to a vacuum polarization. Unfortunately, this state is unstable because 
of the tachyonic mode in the gluon spectrum [21 [3]. This situation is changed at finite 
temperature when the spectrum stabilization happens due to either a gluon magnetic 
mass [6J or a so-called Ao-condensate [5] that is implemented in a stable magnetized 
vacuum. These are the extensions of the Savvidy model to the finite temperature case 
investigated already by the methods of continuum field theory [U [5j [6] . In these ways the 
possibility of the spontaneous generation of strong temperature-dependent and stable 
color magnetic fields of order gH ~ g A T 2 , where g is a gauge coupling constant, T - 
temperature, is realized. The field stabilization is ensured by the temperature and field 
dependent gluon magnetic mass, which serves as a regulator of infrared singularities at 

Such spontaneously created "primordial" color magnetic fields had been generated 
at a GUT scale [7J. They could serve as seed fields responsible for generation of the 
large-scale magnetic fields detected in astrophysical observations [8]. The presence of 
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strong magnetic fields in the early universe is of paramount importance for its evolution. 
In particular, a present day cosmic magnetic field of order B ~ (4 — 5) x 1CT 9 G could 
produce the recently discovered Wilkinson Microscopic Anisotropy Probe (WMAP) 
anomaly [9]. 

In Refs. j6j [TO] the spontaneous creation of the chromomagnetic fields was observed 
in SU(2)- and S'f/(3)-gluodynamics within the one-loop plus daisy resummations. In 
Ref. [11] the chromomagnetic condensate of the same order was calculated in a stochastic 
QCD vacuum model by comparison with some data of lattice simulations. In Ref. 
[T2] the spontaneous generation of the chromomagnetic field at high temperature was 
investigated in a lattice formulation of S'f/(2)-gluodynamics. In Refs. [13] the response 
of the vacuum to the influence of strong external fields at different temperatures was 
investigated and it has been shown that confinement is restored when the strength of 
the external field is increased. 

Other condensate, which formation at T ^ was investigated in continuum field 
theory [H], is the zero component of the gluon electrostatic potential A = const. It 
acts to stabilize the QCD vacuum magnetized state, as it was discussed in Ref. [5]. In 
this paper, however, the magnetic field was considered as external one. Whether or not 
the actual value of the A generated in the deconfining phase is sufficient to remove the 
tachyonic instability was not estimated. 

The main goal of present paper is to investigate the common generation of the 
chromomagnetic field and Ao-condensate in lattice simulations. To incorporate the 
chromomagnetic field on a lattice we use the method developed in [12j. Instead of the 
field strength, which is quantized, the magnetic fluxes are considered as the objects 
to be investigated. The flux takes continuous values. Therefore the minimization of 
free energy of the flux can be done in a usual known procedure. The values of the 
strength of spontaneously generated magnetic field and the values of the Polyakov loop 
for corresponding states were obtained from Monte Carlo (MC) simulations. Then the 
value of A was derived by using a procedure developed in Ref. [15] . 

The paper is organized as follows. In section 2, necessary information about the 
chromomagnetic fluxes on a lattice is educed and the results of calculations are given. 
In section 3, the investigation of the effect of combined generation of chromomagnetic 
field and Ao-condensate is presented. The final section is devoted to discussion. 

2. Chromomagnetic fields on a lattice 

In continuum, to determine the spontaneous generation of a magnetic field one has to 
minimize the effective potential (or free energy) of the background field. The background 
field is introduced by splitting of the gauge field potential A^ into the quantum A^ and 
classical A^ parts: A^ = A^ + A^. We choice the potential A^ a = (0, 0, Hx 1 , 0)8 a3 
that corresponds to a constant chromomagnetic field directed along the third axis in the 
Euclidean and color space. 

We relate the free energy density of the flux to the effective action according to the 
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Figure 1. The plaquette presentation of the twisted boundary conditions. 

definition, 

F(<p) = §{<p) - 5(0), (1) 

where S(<p) and 5(0) are the effective lattice actions with and without chromomagnetic 
field, correspondingly, tp is the field flux. 

To detect the spontaneous creation of the field it is necessary to show that free 
energy has a global minimum at non-zero flux, ip m i n ^ 0. 

In what follows, we use the hypercubic lattice L t x [L t < L s ) with the hypertorus 
geometry; L t and L s are the temporal and the spatial sizes of the lattice, respectively. 
In the limit of L s — > oo the temporal size L t is related to physical temperature. The 
standard Wilson action of the 577(2) lattice gauge theory can be written as 

SV = /?EE [l-ScU*)]; (2) 

Uf, u (x) = U ll (x)U u (x + afi)Ul(x + au)Ul(x), (3) 

where (3 = 4/g 2 is the lattice coupling constant, g is the bare coupling, U^(x) is the link 
variable located on the link leaving the lattice site x in the /i direction, U^x) is the 
ordered product of the link variables. The effective action 5 in (JT]) is the Wilson action 
Sw averaged over the Boltzmann configurations produced in the MC simulations. 

The lattice variable U^x) can be decomposed in terms of the unity, J, and Pauli, 
o~j, matrices in the color space, 

U fl (x)=IUZ(x)+ia j Ul(x). (4) 

The chromomagnetic flux ip through the whole lattice was introduced by applying 
the twisted boundary conditions [16j (Fig. 1). In this approach the edge links in all 
directions are identified as usual periodic boundary conditions except for the links in 
the second spatial direction, for which the additional phase (p is added. The magnetic 
flux (p is measured in angular units and can take a value from to 2tc. 

The twisted boundary conditions (t.b.c.) for the components of (j4j) are 

o j £7° cos <p - C/3 s in <p for x = {x Q , x u L s , x 3 } and (i = 2, 

[f/° for other links, [ ) 
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TI 3( \ j sin <p + cos y? for x = {x Q ,xi, L s ,x 3 } and \i = 2, 
U vW"\ Ul for other links. K ) 

The relations (j5J) and (jHJ) have been implemented into the kernel of the MC 
procedure in order to produce the configurations with the chromomagnetic flux ip. Thus, 
the flux ip is taken into account while obtaining a Boltzmann ensemble at each MC 
iteration. 

The MC simulations are carried out by means of the heat bath method. The 
lattices 2 x 8 3 , 2 x 16 3 and 4 x 8 3 at (3 = 3.0, 5.0 are considered. These values of 
the coupling constant correspond to the deconfinement phase and perturbative regime. 
To thermalize the system, 200-500 iterations are fulfilled. At each working iteration, 
the plaquette value (j3J) is averaged over the whole lattice leading to the Wilson action 
((2j). Then the effective action is calculated by averaging over the 1000-5000 working 
iterations. By setting a set of chromomagnetic fluxes ip in the MC simulations we 
obtain the corresponding set of values of the effective action. The value of the condensed 
chromomagnetic flux <p m in is obtained as the result of minimization of the free energy 
density <^ in the ip. 

The spontaneous generation of chromomagnetic field is the effect of order ~ g 4 
[6]. The results of MC simulations show the comparably large dispersion. So, the large 
amount of the MC data is collected and the standard x 2 - m ethod for the analysis of data 
is applied to determine the effect. We consider the results of the MC simulations as 
observed "experimental data". 

The effective action depends smoothly on the flux ip in the region ip ~ 0. So, the 
free energy density can be fitted by a quadratic function of ip, 

F((p) = F min + b(ip - ipmin) 2 - (7) 

In Eq.(J3), there are three unknown parameters, F min , b and ip m i n . ip m in denotes 
the minimum position of free energy, whereas the F min and b are the free energy density 
at the minimum and the curvature of the free energy function, correspondingly. 

The value ip m in is obtained as a result of minimization of the x 2 -function 

V 2 (F h m \ _ \- (Fmin + K<Pi ~ ip min ) 2 ~ F(ipi)) 2 
X \ r mini 'J i Ymin) 2—j F){F{ip^) ' 

D{F{ipi)) = £ (F( ^~_^ m)2 , (9) 

where ipi is the array of the set of fluxes and D(F(ipi)) is the data dispersion, which can 
be obtained by collecting the data into bins (as a function of <p)\ n bin is a number of 
points and and F^n is a mean value of free energy density in the considered bin. As it 
was determined in the data analysis, the dispersion is independent of ip. The deviation 
of ipmin from zero indicates the presence of the spontaneously generated field. 

The fit results are given in the Table 1. As one can see, ip m in demonstrates the 
2<r-deviation from zero. 

The 95% C.L. domain of parameters F min {b for the right figure) and ip m in is 
represented in Fig. 2. The black cross marks the position of the maximum-likelihood 
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Figure 2. The 95% C.L. area for the parameters F m i n and ip m in, determining the free 
energy density dependence on the flux <fi m in on lattice 2 x 8 3 for j3 = 3.0 (left figure). 
The 95% C.L. area for the parameters F m i n and b , determining the free energy density 
dependence on the flux <p m i n on lattice 2 x 16 3 for f3 — 3.0 (right figure). 



values of F min {b for the right figure) and (p m in- It can be seen that the flux is positively 
determined. The 95% C.L. area becomes more symmetric with the center at the F m i n , 
b and (p m in when the statistics is increased. This also confirms the results of the carried 
out fitting. 



3. The spontaneous vacuum magnetization and Ao-condensate 

In this section we investigate an effective potential, taking into consideration the effects 
of non-trivial A -condensate and chromomagnetic field |15J . This is realized by means 
of calculation of the partition function in a general covariant background field gauge 
with the background field providing both the chromomagnetic field and the non-trivial 
Polyakov loop. Then the data obtained in MC simulations are substituted in this 
effective potential. If the minimum value of it is negative and the potential is real, 
the common generation of condensates happens. 

Within the used imaginary time formalism, the temporal direction in the Eucleadian 
space is periodic, with period (3 = 1/T, and the Polyakov loop Pi, defined as the time- 
ordered exponential, 



Pl(x) = Texp 



ig <j> drA (x, r) 



(10) 



is a proper order parameter to describe the deconfinement phase transition. It is specified 



Table 1. The values of the generated fluxes (p m in for different lattices. 
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by a constant A field, given in the fundamental representation by 



A n = ^ 



11) 



90 2' 

where < < 2ir, r 3 is Pauli matrix. 

The trace of the Polyakov loop in the fundamental representation is 

Tr F (P L ) =2cos(0/2). (12) 

Unless the Z(2)symmetry is spontaneously broken, the variable has to have the value 
7r corresponding to Tr^(P^) = 0. 

The one-loop contribution to the free energy has the form [TJ] 



yd) = vi l) + vf> 



(13) 



E 1 



EE 

m=0 n,± 



d 3 k 
lgH 



log 



/3 2tt 7 2tt 



° ^ ^ /5 2vr 



m=0 n,± 



rffc 3 
Z7T 



(o; n -p 2 + 2^F(m+-±l) 
K + |) 2 + 2^(m+i±l) 



/c 2 



k\ 



+ 



where uj n = 2im/f3 is the Matsubara frequency, and the sum over n is over all integer 
values. The sum over m is the sum over the Landau levels. 

By applying a standard technique of Schwinger to express the logarithms of ( [TBI 
as an integral and performing the high-temperature expansion of the effective potential, 
we obtain 

(gH) 2 (gH) 3 / 2 ™2 2k B 2 , 



E 

i 



^" r k=2 



P 2 gH 



(2k)\ 

2fc-3/2 / nI j\2 



Y{2k - 3/2) x 



(14) 



(3 2 gH+ (0-2vr/) 2 



(gH) 2 , x 2tt 2 
Ky ' (3 - 4 7 ) - 



+ 



3/? 4 ' 12vr 2 /? 4 

(ffif/3 2 + 2 ) 3 / 2 
3^ 
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27T/3 2 

S 2 # 2 log 



gH(3 2 



24tt 2 
+ 0^ + 
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S 2 # 2 



12%^gHp 2 + 2 



4lT 



An 2 
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9 2 H 2 log 



Lit 



12tt 2 



^ 3 /3 2 C 3,-£ ^ 3 /3 2 C 3,£ 



2tt 



3/92/ 



2- 



+ 0(/5 2 (^) 3 ). 



24tt 4 192vr 4 192vr 4 

Here C\ = —0.01646, 7 ~ 0.577216 is Euller's constant, ((s) is the Riemann Zeta- 
function, B n are the Bernoulli numbers. 

We performed the described above lattice calculations and have obtained the values 
for the spontaneously generated field strength and the Polyakov loop of corresponding 
states. Then we considered the expression (|14p as a function of these parameters 
and studied its minimum. We have determined that for the temperatures 1/(3 = 
100 - 500GW, the field strengths H = 10 22 - 10 24 G, and = 2.0-3.11 the one- 
loop effective potential (fill) is negative, V^(H, 0, (3) < 0. Hence it follows, either 
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gauge field or A condensate are spontaneously generated. Moreover, in the minimum 
of the effective potential for used values of the parameters the ratio ImV^ / ReV^ is 
of order ~ 10~ 10 , that is practically zero for numeric calculations. We conclude that the 
common generation of the constant Abelian chromomagnetic field and the non-trivial 
Polyakov loop results in the stable magnetized vacuum state. Of course, the one-loop 
potential does not exhibit the complete effective potential. However, it gives the main 
contribution which contains a larger imaginary part. If one adds the second next-to- 
leading contribution coming from ring diagrams, this part is cancelled [6]. So, the check 
of stabilization at one-loop level is important. 

4. Discussion 

We have investigated the effect on common spontaneous generation of chromomagnetic 
field and A -condensate in SI/(2)-gluodynamics at high temperature. The first step 
was to show the possibility of the spontaneous generation of chromomagnetic field at 
high temperature [12J . The obtained results are in a good agreement with that of 
derived already in the continuum quantum field theory 0, [17] and in the lattice data 
analysis pXj . The actual values of the chrmomomagnetic field strength and the Polyakov 
loop were obtained from the same MC simulations. They allow to conclude that these 
condensates, chromomagnetic field and A , have to be present in the deconfinement 
phase of QCD. 

The developed approach joins the calculation of free energy functional F((f>) and the 
consequent statistical analysis of its minimum positions for various temperatures and 
flux values. In this way the spontaneous creation of condensates is realized in lattice 
simulations. 

Other important result of the present work is that the spontaneously magnetized 
state is stable. This indicates that a true vacuum at high temperature is formed from 
these condensates. 

Since the field configuration with constant fields is gauge non-invariant, one could 
consider this configuration as a domain. A complete structure of whole space can be 
derived by using the requirement of gauge invariance for the space including domains 
with different orientation. This problem we left for the future. 

We would like to conclude that in the deconfinement phase the condensates have to 
influence various processes that should be taken into consideration to have an adequate 
concept about this state. 
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